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OBSERVACIONES AL ALUMNO 

Al comienzo de cada texto encontrará lo siguiente. 

1. UNA PRESENTACION DE NUESTROS OBJETIVOS en la unidad, la cual puede usar para comprobar si ha completado el 
curso satisfactoriamente. 

2. UN DIAGRAMA ESTRUCTURAL que le muestra cómo están relacionadas las diferentes secciones de la unidad. El diagrama 
consta de: 
a) cuadros rojos que muestran la secuencia principal. 
b) cuadros negros punteados (a la izquierda del rojo) que le muestran lo que esperamos ya sepa antes de fijarse un estudio 

detallado. 

e) cuadros negros (a la derecha del rojo) que le ayudan a escoger temas que puede posponer u omitir si tiene »0co tiempo. 

3. UN GLOSARIO de las definiciones, notaciones, etc., usadas en la unidad como referencia inmediata. 

4. UNA BIBLIOGRAFIA de los libros que puede consultar cuando haya acabado el trabajo de la unidad. 

5. CONCLUSIONES Y RESUMEN del trabajo realizado al final de cada texto. 


NUMERACION E INDICADORES 

Dentro de cualquier sección nos referimos a los temas por medio de un número. Al referirnos al mismo tema en otra sección, citamos 
cuatro números: número de la unidad, número del capítulo o parte, número de sección y número del tema. Por ejemplo, “Ecuación 3” 
si aparece en la misma sección, y “Ecuación 4.3.1.3”, si está en otra sección. Los tres primeros números se citan en la parte superior 
de la página. Creemos que este es un sistema claro y fácil para consultar cualquier tema. Los indicadores serán una guía valiosa 
y le ahorrarán tiempo. 

En la margen derecha del texto hemos tratado de indicar la clase del material. Los términos usados son tan claros que no necesitan 
explicación. Por ejemplo, “Discusión” significa que una vez que haya entendido el tema, no necesita de nuevo recurrir a esta parte 
del texto. Estos INDICADORES tienen por objeto ayudarlo a localizar rápidamente las partes importantes cuando esté repasando 
o tenga que omitir parte del texto por falta de tiempo. 

También hemos tratado de indicar la importancia de los temas por el sistema de estrellas. 

*** Esto es muy importante. Usted debe entenderlo claramente. En algunos casos lo usará con tanta frecuencia, que se familiari- 
zará con ello sin necesidad de memorizarlo deliberadamente. El punto principal de este curso es no aprenderse los datos de me- 
moria; nosotros creemos que es más importante que entienda el trabajo y sea capaz de aplicarlo en la solución de los problemas. 
** Este es de menor importancia, pero tiene material al cual nos referiremos más adelante. 

* Esto generalmente no lo guía a trabajos adicionales en este curso y lo puede omitir si tiene poco tiempo. 

M Indica el final de cada ejemplo, ejercicio o solución. 
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Objetivos 


El principal objetivo de esta unidad es discutir el problema de la 
solución de un sistema de ecuaciones lineales simultáneas y, en par- 
ticular, discutir la exactitud de los resultados y la eficiencia del 
método usado. 


Después de haber trabajado esta unidad, debe estar en capacidad de: 


(i) explicar el significado de los siguientes términos: 
vector error, 
sistema de ecuaciones de mal comportamiento, 
matriz casi singular, 
norma de un vector; 


(ii) determinar la eficiencia del método de eliminación de Gauss o 
un método de matriz inversa, determinando el número de opera- 
ciones aritméticas requeridas; 


(i1i) reordenar una ecuación matricial dada 
Ax =b 


de manera que se garantice la convergencia de un método de so- 
lución iterativo; 


(iv) determinar si un sistema (sencillo) dado de ecuaciones simul- 
táneas es de mal comportamiento. 


iv 


Diagrama estructural 
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Una solución formal 
de un sistema de 
ecuaciones, 

28.1.1 


Método de 
eliminación de 
Gauss, 

28.1.2 


Método para 
invertir una 


matriz, 
28.1.3 


Método iterativo, 


28.2 


Sistemas de 
ecuaciones de mal 


comportamiento, 
28.3 
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Ejercicio opcional 


de programación, 
28.2.2 


Glosario 


MATRIZ CASI 
SINGULAR 


MATRIZ DE LOS 
COEFICIENTES 


MATRIZ RALA 


METODO DIRECTO 


METODO INDIRECTO 


NORMA DE UN 
VECTOR 


SISTEMA DE 
ECUACIONES DE MAL 
COMPORTAMIENTO 


VECTOR ERROR 


vi 


Una matriz casi singular es una matriz no sin- 
gular tal que pequeñas variaciones en los ele- 
mentos producen una matriz singular. 


La matriz de los coeficientes del sistema de 
ecuaciones 


011X1 + 412X2 +: +41,X, =b, 


071X1 + 022X2 + *** + Q2pX, =b, 


Am1X1 + Qm2X2 + *** + QAmpXn = Bn 


es 
411 412 Ain 
421 422 02n 
Ami Am2 Arm 


Una matriz rala es una matriz donde los ele- 
mentos predominantes (por la frecuencia con 
que aparecen) son ceros. 


Un método directo para resolver un sistema de 
ecuaciones lineales simultáneas es un método 
que produce un resultado explícito de la so- 
lución. 


Un método indirecto es, en esta unidad, un mé- 
todo iterativo. 


Una norma de un vector es una medida numé- 
rica del “tamaño” del vector. 


Un sistema de ecuaciones de mal comporta- 
miento es un sistema donde ocurre que peque- 
ñas variaciones en los valores de los coeficien- 
tes producen grandes cambios en los valores de 
los elementos del conjunto solución. 


El vector error de un sistema de ecuaciones 
lineales es 


(vector solución 'estimado) — (vector solución) 
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Notación 


Los símbolos están dados en el mismo orden en que aparecen en el 
texto. 


Á Una matriz de coeficientes no singular. 
Xx Un vector. 
E El conjunto de todos los números complejos. 
AOS La inversa de A. 
Xi Un elemento del vector x. 
0; Un elemento de la matriz de coeficientes A. 
x0 La r-ésima aproximación del vector solución en un método 
iterativo. 
x” Un elemento de x!”. 
X El vector solución exacto. 
gr El vector error: 
(5% — X) 
en Un elemento de el”. 
lal Una norma del vector a definida por 


llall = la] + laz] + --- + la), 
donde a es un vector n X 1 de elementos a;, i =1,...,n. 


i.e. Una abreviatura para significar “esto es” o “es decir”. 
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L. Fox, An Introduction to Numerical Linear Algebra (Claredon Press, 
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Este libro es un estudio muy amplio de la teoría y la práctica de los 
métodos directos e indirectos de resolver ecuaciones lineales simul- 
táneas y de análisis de error asociado. Va más allá de lo que se pro- 
pone este libro, pero se hace referencia a varios puntos para beneficio 
del estudiante que desea avanzar en el tema. 


A. Ralston, A First Course in Numerical Analysis (McGraw-Hill, 1965). 


Este libro estudia muchos aspectos del análisis numérico así como los 
que se incluyen con la solución de las ecuaciones lineales simultáneas. 
El Capítulo 9 es una buena lectura como estudio adicional sobre solu- 
ciones de ecuaciones lineales simultáneas, especialmente desde el 
punto de vista del análisis de error. 


Prepárese: Usted encontrará ambos libros muy difíciles si sus conoci- 


mientos en este tema se basan únicamente en el material de esta 
unidad. 
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28.0 Introducción 


Esta unidad considera detalladamente un aspecto del trabajo discu- 
tido en la Unidad 26, Algebra lineal HI. Vamos a analizar algunos 
métodos prácticos para resolver un conjunto de ecuaciones lineales 
simultáneas, y algunas de las dificultades que se presentan. Exami- 
naremos este problema como un vehículo para la discusión de algunos 
de los métodos y prácticas del cálculo numérico: es otro paso en el 
aspecto del análisis numérico de este curso, que comenzó en la Unidad 
2, Errores y exactitud. 


Comenzaremos por establecer el problema con toda precisión. 


La ecuación matricial 
Ax = D, 


donde x y b son vectores columna de n elementos y A es una matriz 
n X n, representa un sistema de n ecuaciones en n variables (in- 
cógnitas). Si A es una matriz no singular, esto es, si su rango es n, 
entonces la ecuación matricial tiene una solución única. Este es el 
único tipo de ecuaciones que consideraremos en esta unidad porque 
uno de nuestros objetivos es examinar cómo se llega a la solución, 
particularmente en el caso de las matrices grandes. Así, pues, supon- 
dremos que existe la solución única; nuestro problema es determi- 
narla. 


Si ha resuelto ecuaciones simultáneas anteriormente, probablemente 
ha resuelto sistemas de 2 ó 3, o quizá 4, ecuaciones. En su trabajo, se 
esforzó por evitar errores y, como comprobación final, remplazó sus 
respuestas en las ecuaciones originales. A usted no le interesaba com- 
parar unos métodos con otros para ver si alguno era más rápido que 
los demás; solamente se preocupaba por la exactitud de la respuesta, 
en el sentido de que éstas verificaran el sistema original. Cuando se 
tienen sistemas de varios miles de ecuaciones, el tiempo puede ser un 
factor importante y veremos que hay más exactitud que cálculos arit- 
méticos correctos. Cuando resuelve sus ecuaciones, es casi seguro 
que utiliza métodos directos, esto es, un método que, paso por paso, lo 
lleva a la respuesta, la cual aparece en el último paso. Por ejemplo, un 
método directo es el de eliminación de Gauss que analizamos en la 
Unidad 26. Lo examinaremos nuevamente, con otros dos métodos di- 
rectos, en la Sección 28.1. 


Cuando resuelve un sistema de ecuaciones pequeño, probablemente 
no pienza (ni ha pensado nunca) en un método iterativo (es decir, en 
un método en el cual vamos obteniendo, en cada paso, una aproxima- 
ción de la solución); la razón es que el método directo es rápido y efec- 
tivo. En los sistemas grandes, con ciertas características, pueden ser 
útiles los métodos iterativos. Examinaremos estos métodos en la 
Sección 28.2 y veremos cómo se relacionan con las ideas de los méto- 
dos iterativos de solución de ecuaciones (descritos en la Unidad 2, 
Errores y exactitud) y con algunos de los conceptos sobre espacios 
vectoriales vistos en las unidades anteriores. 


Por último, llegamos al problema de la exactitud. Frecuentemente, los 
datos que utilizamos para construir las ecuaciones no son exactos y, 
en consecuencia, los resultados eventuales pueden contener errores, 
no solamente por causa del redondeo de los errores que vamos hacien- 
do durante el proceso, sino por las inexactitudes que se van propa- 
gando desde el inicio. En ciertos casos especiales, esta inexactitud 
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hace los resultados casi inaceptables. Es este aspecto de la solución 
de los sistemas de ecuaciones, llamados de mal comportamiento, el 
que estudiaremos en la Sección 28.3. 


28.1 METODOS DIRECTOS 


28.1.0 Introducción 


A los matemáticos les gusta poder registrar una solución explícita de 
una ecuación. La solución, tantas veces citada de la ecuación cua- 
drática 


ax?+bx+c=0, 


a saber 


_—b+ yb? -— 4ac 


> 2a j 
es un ejemplo de esto. La variable deseada está a la izquierda de un 
signo “igual a”, y las letras a la derecha del signo “igual a” se pueden 
sustituir por números tomados de un ejemplo específico dado. La so- 
lución consiste en los dos números que están aplicados al 0 en la 
función 


x—ax?+bx+c  (xeC). 


Del mismo modo, podemos escribir una solución explícita de la ecua- 
ción matricial 


Áx =b, 
donde A es una matriz no singular, en la forma 
x=A"'!b, 


que representa el vector solución, que es el único vector aplicado a b 
en la aplicación 


AA (x e R”). 


Cualquier método que utilice este tipo de fórmula directa para obtener 
la respuesta se llama método directo. No obstante, no todos los méto- 
dos directos utilizan fórmulas para calcular la respuesta. Por ejemplo, 
el método de eliminación de Gauss (descrito en la Unidad 26) es un 
método directo pero, sin embargo, no utilizamos en él una fórmula 
para hallar la respuesta, sino un procedimiento que, paso por paso, va 
desde los datos hasta la respuesta (estos pasos se podrían especificar 
mediante un algoritmo). Esta es la característica de los métodos direc- 
tos (en oposición a los métodos iterativos, que discutiremos más 
tarde): es un método para obtener la respuesta de un problema, a 
partir de unos datos, mediante un procedimiento que va paso por 
paso. La respuesta (o una aproximación de la respuesta) se obtiene 
únicamente al concluir todo el proceso; no se obtienen aproximaciones 
de la respuesta durante el desarrollo del proceso. 


Para calcular el vector solución de un sistema de ecuaciones, median- 
te la fórmula x = A”'b, tenemos que determinar A7! y multipli- 
carla, a la derecha, por b. En la Sección 28.1.3 examinamos la eficien- 
cia de este proceso particular. También compararemos su eficiencia 
con la eficiencia de otros dos métodos: el método de eliminación 
de Gauss, que examinaremos en la Sección 28.1.2, y otro método direc- 
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* * 


Definición 1 
hk Rx * 


to, que veremos en la Sección 28.1.1. Por último, debemos destacar 
que los tres métodos de solución son equivalentes puesto que la solu- 
ción es única. Esto significa que, si los escribiéramos explícitamente 
mediante fórmulas, podríamos derivar una fórmula de otra cualquiera, 
por medio de convenientes manipulaciones algebraicas. En esta uni- 
dad nos interesaremos en comparar unos con otros los métodos opera- 
cionales y los tiempos que requiere cada uno. 


28.1.1 Un método explícito 


Mediante operaciones adecuadas (por ejemplo, mediante la elimina- 
ción de x2 en las ecuaciones y el consiguiente despejo de x,), po- 
demos escribir la solución de 


411X1 + 412X> =b, 
471X1 + 477X, = b, 
en la forma 


_ b,a3, — b3a,) oe 411b3 — a31b, 


X1 2 


= , E , 
411422 — 412421 411422 — 412421 


supuesto que 
4,1472 — 41247, H0. 


(Supondremos que todos los coeficientes a, aja, ... y bi, ba, 
utilizados aquí son números reales.) En este caso, la solución encon- 
trada es el conjunto de los elementos componentes del vector x (único) 
que está aplicado a b en la aplicación 


X>Ax  (xER?) 
Así, en un paso, suponiendo que se verifica la condición, hemos ex- 


presado formalmente la solución de todos los pares de ecuaciones 
simultáneas en dos variables, que tienen solución única. 


Podemos proceder a buscar este tipo de solución para un sistema de 
tres ecuaciones lineales simultáneas: 


411X1 + 412X) + 013X3 =b, 
471X, + 077X2 + 473X3 = b, 
431X1 + 437X2 + 433X3 = b, 
cuyo conjunto solución se puede expresar como 


A b,(472433 — 473432) — by(4,2433 — 4,3432) + ba(a,2473 — 4,3422) 
4, ¡(4,2433 — 473432) — 421(4,2433 — 4,3432) + a3,(a, 7473 — 4,3472) 


X1 


conjuntamente con las expresiones similares para x2 y x3 (que 
tienen el mismo denominador), supuesto, desde luego, que el denomi- 
nador no es igual a cero. (Si tiene tiempo suficiente, quizá querría 
verificar por sí mismo estas expresiones, resolviendo el sistema de 
ecuaciones. Sin embargo, no deje que las excesivas manipulaciones 
algebraicas lo demoren demasiado.) 


También podemos ver aquí que, por una simple sustitución de letras 
por números, es posible resolver cualquier sistema de la forma anterior 
que tenga solución única. Podríamos, en teoría, resolver un sistema 
de ecuaciones lineales simultáneas de cualquier orden, desarrollando 


3 


MB 28.1.0, 28.1.1 


28.1.1 


Tema principal 
VW 


las fórmulas apropiadas. ¿Por qué, pues, no nos detenemos aquí y nos 
limitamos a resolver los sistemas de ecuaciones simultáneas utili- 
zando fórmulas convenientes? La razón es que el gasto de tiempo y 
esfuerzo es sumamente grande. 


Consideremos el número de multiplicaciones y divisiones comprendi- 
das en la solución de tres ecuaciones lineales simultáneas. Nos referi- 
mos a las multiplicaciones y divisiones porque estas operaciones 
requieren más tiempo que las adiciones y sustracciones, ya sea que 
se efectúen a mano o por medio de alguna máquina. 


Las expresiones que corresponden a x¡, x2 y x3 tienen el mismo 
denominador, de modo que basta calcularlo una sola vez. Cada nume- 
rador tiene la misma forma que el denominador, y hay 3 numeradores. 
Así, necesitamos calcular 4 expresiones de la forma 


4, ¡(4,2433 — 4,343>,) — 471 (4,7433 — 4,3432) + a3,(4,,473 — 4,3422). 


Cada paréntesis contiene 2 productos. Por consiguiente, la evaluación 
de toda la expresión requiere 9 multiplicaciones. 


El número total de multiplicaciones es 4 X 9 = 36 
Número de divisiones = 1) 
Por tanto, el número de operaciones, que llevan tiempo, es = 39 


Usted debe haber observado que los términos encerrados en los parén- 
tesis en cualquiera de los mumeradores (en el de x, en nuestro 
ejemplo) son iguales a los términos que aparecen dentro de paréntesis 
en el denominador. Aprovechando ese hecho, podemos ahorrarnos 6 
multiplicaciones y tener solamente 33 operaciones, que llevan tiempo. 


Ejercicio 1 


¿Cuál es el número total de operaciones de multiplicar y dividir que se 
requieren para resolver un sistema de dos ecuaciones lineales simul- 
táneas mediante las fórmulas dadas en el texto. al 


Se puede deducir una expresión general del número de operaciones 
de multiplicar y dividir que se requieren para la solución de un siste- 
ma de n ecuaciones con n incógnitas, por este método; viene dada 
por la sucesión 


u,=8 


e n> 2. 


Se puede demostrar que, para valores grandes de n, el n-ésimo térmi- 
no de esta sucesión, u,, es aproximadamente igual a 


Liam <a 


Esto nos permite construir la tabla que damos a continuación. La 
última columna de la tabla nos da una ligera idea del tiempo que ne- 
cesitaría emplear un computador automático para resolver el proble- 
ma si, por ejemplo, cada operación requiere un microsegundo (1 micro- 
segundo = 107”* seg; se denota por seg). 


Esto muestra el tremendo crecimiento del número de operaciones y 
del consiguiente tiempo requerido, cuando n crece. Se necesita, ob- 
viamente, un método más eficiente. 
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Número de Número de Tiempo 
ecuaciones simultáneas| operaciones requerido 


8 
33 
168 
=6,2 x 10” 
=1,6x 101% 
=6,9 x 102570 


Sección opcional sobre determinantes Material opcional 
* 


Si desea, puede omitir esta sección, que no es esencial para el des- 
arrollo de esta unidad ni para el desarrollo del curso. La incluimos 
porque ella analiza brevemente los determinantes, que, en un tiempo, 
estuvieron íntimamente ligados a la solución de las ecuaciones linea- 
les simultáneas. 


La expresión 
411422 — 412421 
se llama determinante de la matriz 
a11 A 
A) = , 
421 422 
y se expresa por 
A11 412 
421 422 


o det Az o |42|, de modo que la solución de dos ecuaciones en dos 
variables se puede escribir como 


x, = ——— etc. 


- Del mismo modo, la expresión 
4, ¡(47,433 — 423432) — 4,1 (4,2433 — 4,3432) + a31(4,,473 — 4,3422) 


es el determinante de la matriz 


411 42 413 
A3=|4,, 47 43], 
431 432 433 


y se expresa por 


411 412 413 
421 422 43 


a ¿8 08 
431 432 433 (continúa en la página 6) 
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Solución 1 Solución 1 


8 5 


(viene de la página 5) 


o det Az o |A3|, de modo que la solución de tres ecuaciones con 
tres incógnitas se puede escribir como 


b, a, 413 
bz 42, 423 


b3 432 433 
411 412 413 


etc. 


El determinante de una matriz cuadrada es un cierto número que está 
asociado a la matriz. Esencialmente, su uso, en este nivel matemáti- 
co, se reduce a ser una abreviatura para expresar las soluciones de 
ecuaciones lineales simultáneas. Puesto que el uso real de las fórmu- 
las para encontrar las soluciones (fórmulas que son equivalentes a la 
evaluación mediante determinantes), se hace muy complicado cuan- 
do el sistema consta de más de 3 ó 4 ecuaciones lineales simultáneas, 
no hemos querido que los determinantes sean una parte integral del 
curso. El número de ecuaciones óptimo, por así decirlo, para el uso de 
los determinantes es solamente tres. 


28.1.2 El método de eliminación de Gauss 28.1.2 
En la Unidad 26 introdujimos el método de eliminación de Gauss para Tema principal 
la solución de los sistemas de ecuaciones. El objetivo del método es *.* 


producir un conjunto equivalente de ecuaciones que son más fáciles 
de resolver que las del conjunto original. Considere ahora este método 
y compruebe su eficiencia en términos del número de operaciones que 
necesita para llegar a la respuesta, puesto que, para sistemas grandes, 
el método discutido en la Sección 28.1.1 es, evidentemente, inacep- 
table. (Si alguien, a principios de este siglo, hubiera contado con un 
computador automático y se hubiera dispuesto a resolver por ese mé- 
todo un sistema de 100 ecuaciones simultáneas, no habría llegado 
todavía a hacer mella en el problema: aún no tendría realizado ni si- 
quiera el 107 !%%% de los cómputos necesarios.) 


Ejercicio 1 Ejercicio 1 
sE , E Ñ , (3 minutos) 
Resuelva el siguiente sistema de tres ecuaciones por medio del méto- 


do de eliminación de Gauss con sustituciones regresivas, en el orden 
mostrado. Realice las operaciones en un papel y, después, utilice los 
espacios destinados a los distintos pasos del proceso. Los números 
deben expresarse en forma decimal, y no como fracciones. Escriba (en 
la columna de la derecha) el número de divisiones y de multiplica- 
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ciones que ocurren. ¿Cuál es el número total de estas operaciones 
que se requieren para obtener la respuesta correcta? 


Sx1 + 2x, > 2x; = 8 (1) 
3x1 + 6x, > 4x3 5 1 (2) 
2x1 + 4x, el 2x, = 1 (3) 


Número de multiplicaciones 
y/o divisiones 


Paso en el proceso 


Elimine x, de las ecuaciones (2) y 


El factor de la ecuación 
(2) es 


El factor de la ecuación 
(3) es 
Sx 1 2x, + (-2) 


pu 
= 


> 
w 
Il 


00-90 


> 
w 
| 


(4) 


a 
[91] 
= 

— 

un 

= 


Ahora elimine x2 de la ecuación 


El factor de la ecuación 
(5) es 


Sx, + 2x2 + (-2) X3 = (1) 
(4) 
3 = (6) 


Por consiguiente, 


e 
w 
Il 


5 
mi 
Ju) 


y la sustitución regresiva en la ecuación 


(4) da 


X2 = 


y de la ecuación (1) 


ON 


Número total de divisiones y multiplica- 
ciones 


ps (continúa en la página 8) 
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Solución 1 Solución 1 


Número de multiplicaciones 
Paso en el proceso y/o divisiones 


Elimine x, de las ecuaciones (2) y (3). 


El factor de la ecuación (2) es —0,6 |1 división 


El factor de la ecuación (3) es —0,4 |1 división 
Sx, + 2x2 + (-2) X3 = 8 (1) 


4,8 x+ -28  x; —3,8 (4) |3 multiplicaciones 


3,2 Xx2+ —112  X3 —2,2 (S) |3 multiplicaciones 
Ahora elimine x2 de la ecuación (5). 

El factor de la ecuación (5) es —0,66 |1 división 

Sx, + 2x2 + (-2) X3 


4,8 X2 + —2,8 X3 


0,66 x3= 2 multiplicaciones 


Por consiguiente, 


X3 = 1 división 


y la sustitución regresiva en la ecuación 
(4) da 
x= —0,5 1 multiplicación y 1 división 


y de la ecuación (1) 


=> 2 2 multiplicaciones y 1 división 


Número total de divisiones y multiplica- 
ciones 17 


Observe que el número de multiplicaciones y divisiones es, aproxima- 
damente, la mitad del número de operaciones que requiere el método 
considerado en la Sección 28.1.1. E 


(viene de la página 7) 


Ahora trataremos de determinar el número de operaciones que se Tema principal 
requieren para resolver un sistema general de ecuaciones por el mé- *o* 

todo de eliminación de Gauss. Hagamos un análisis semejante al que 

hicimos en el último ejercicio. Comparemos, después, el trabajo exigido 

por este método con el trabajo exigido por el de la Sección 28.1.1. 
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Consideremos las dos primeras ecuaciones de un sistema de n ecua- 
ciones: 


4,¡X1 + 419X2 + +** + 01, =b, 
471X1 + 077X2 + :** + Oa,X, =D) 


Supondremos que a, es diferente de O. Si no fuera diferente de O, 
tendríamos que alterar el orden de las ecuaciones hasta encontrar un 
coeficiente de x, que sea diferente de 0. Cuando las operaciones se 
efectúan a mano, esta reordenación de las ecuaciones no implica un 
cambio apreciable del tiempo requerido, pero en los cálculos auto- 
máticos tendrá, por supuesto, que tenerse en cuenta. Para la elimina- 
ción de x, de la segunda ecuación se requiere el cálculo del cociente 
az,/a11 y, después, la formación de los n números 


Así, la formación de la nueva segunda ecuación requiere (n + 1) ope- 
raciones de multiplicación o de división. En el siguiente ejercicio 
tendrá que calcular el número total de multiplicaciones y divisiones 
requeridas. 


Pe 421 a 
Para el próximo paso supondremos que 4», — a a,, es diferente de 0 
11 


porque tendremos que utilizarlo como divisor para obtener el cociente 


: 421 
correspondiente a —=, esto es, 


411 


432 


421 
422 7—Uia 
411 


Si fuera 0, tendríamos que alterar el orden de las ecuaciones hasta 
encontrar un coeficiente de x2 que sea diferente de 0. (¿Qué po- 
dríamos inferir si todos los siguientes coeficientes de x2 fueran 
iguales a 0?) 


En general, supondremos que, en las discusiones que aparecen en 
este texto, es siempre posible efectuar todas las divisiones que van 
ocurriendo. 


Ejercicio 2 Ejercicio 2 
dd Ñ , E (4 minutos) 
Llene la siguiente tabla que está destinada a calcular el número de 


multiplicaciones y divisiones que se requieren para resolver un siste- 
ma de n ecuaciones con n incógnitas por el método de eliminación de 
Gauss. Las fórmulas obtenidas en la Unidad 4, Diferencias finitas, 


a e ED 
r=1 2 
ra Me + Du D 
27 = 6 > 


serán muy útiles. 
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Número de multiplicaciones 
y/o divisiones 
(n — 1)n+1)=n?* —1 


Paso en el proceso 


Elimine x, de todas las ecuaciones 
que siguen a la primera. 


Elimine x2 de todas las ecuaciones 
que siguen a la segunda. 


Elimine x,-, de todas las ecuaciones 
que siguen a la (n — 1)-ésima, esto es, 
de la n-ésima ecuación. 


El total, para las eliminaciones, es 


Ahora efectúe las sustituciones regresivas. 


Para determinar x, de una ecuación tal 
como 


An = Br. 
Para determinar x,- de una ecuación 
tal como 


Yn-1Xn-1 + YnXn = Bn-1 


Para determinar xi 


El total para las sustituciones regresivas 
es 


Por consiguiente, el número total de 
Operaciones es 


Compare su respuesta con la respuesta del ejercicio 1. S 


Ahora estamos en condiciones de poder comparar la solución explícita 
y el método de eliminación de Gauss, desde el punto de vista de la 
eficiencia. Los resultados se dan en la tabla: 


Método de la fórmula Método de eliminación 
El número de de Gauss 
Operaciones es El número de operaciones es 
Número de 1,72n X n! 
ecuaciones (para n grande) 


33 
168 
=6,2x 10” 
=16x 10160 
=6,9 x 102570 


Esta tabla muestra claramente por qué no se utiliza el método de la 
fórmula cuando se trabaja con sistemas grandes de ecuaciones. 


De hecho, en ningún caso el método de la fórmula es más rápido que 
el de eliminación de Gauss. Para establecer una rápida comparación 
de los métodos, cuando n es. grande, diremos que el número de opera- 
ciones que requiere el método de eliminación de Gauss es, aproxima- 


n ; 
3» Porque, cuando n > 100 por ejemplo, el resto de los 
términos afecta al número total de operaciones en menos de un tres 


por ciento. 


damente, 
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Discusión 
* * 
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Solución 2 Solución 2 


Número de multiplicaciones 
Paso en el proceso y/o divisiones y 


Elimine x, de todas las ecuaciones que|  (n— 1)(n + 1) = n—1 
siguen a la primera. 


Elimine xz de todas las ecuaciones que|  (n — 2)n = (n — 1)? — 1 
siguen a la segunda. 


Elimine x,-, de todas las ecuaciones que 
siguen a la (n — 1)-ésima, esto es, de la 
n-ésima ecuación. 


El total, para las eliminaciones, es == 1) 


nín + 1)(2n + 1) 
A 1) Pa 


Ahora efectúe las sustituciones regresivas. 


Para determinar x, de una ecuación tal 
como 

AX = Ba 
Para determinar x,-, de una ecuación tal 


como 
Yn-1%n-1 «+ YnXn = Br=1 


Para determinar x; 


El total para las sustituciones regresivas 
es 


Por consiguiente, el número total de 
Operaciones es 


28.1.3 Un método para calcular la matriz inversa 


En la Unidad 26 describimos un método numérico para encontrar la 
inversa de una matriz n X n; el método es muy semejante al método 
de eliminación de Gauss. Señalaremos que, si tuviéramos que resolver 
varios sistemas de ecuaciones de la forma 


¡Ax =D, 


donde la matriz A es la misma en cada caso, y solamente cambia la 
matriz b, entonces sería de algún valor encontrar A”! y calcular la 
solución mediante la aplicación de la fórmula 


y = A” 1h 
a cada caso. 


Examinaremos ahora el número de operaciones requeridas y veremos 
cuándo este método es de “algún valor”. 


Para recordar el método describiremos los pasos esenciales. Suponga- 
mos que queremos hallar la inversa de la matriz 


2 1 -—1 
6 —1 -—9 
4 3 1 
Escribimos 
2 il =4 | 1 0. ..0 
6-1 -910 10 
> 3 ¡4 0 1 


1 0 0 x x Xx 
0 1 0 x x x 
0 0 1 x x Xx 


donde las cruces representan los elementos de la matriz inversa. He- 
mos omitido la columna de las sumas de verificación: esta columna 
debe aparecer en ambos métodos y, en consecuencia, podemos omitirla 
en una comparación. 


Comenzamos por dividir la primera fila de nuestra disposición de 
números por 2 para obtener un 1 en la esquina superior izquierda. 
Esto implica 3 divisiones. Obtenemos 


E 
| 

6-1 $10 1 0 

e. 3 106 0 1 


Después restamos 6 veces la primera fila de la segunda, y 4 veces la 
primera fila de la tercera. Tenemos 2 X 3 multiplicaciones. Obtenemos 


1 44; 34 0 0 
| 
4 613 1 
E E A E | 
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28.1.3 


Tema principal 
* + + 


Así, para llevar la primera columna a la forma requerida hemos tenido 
que efectuar 9 operaciones, que llevan tiempo. 


En general, para calcular la inversa de una matriz n X n, formamos 
, 

la correspondiente matriz n X 2n “uniendo” a la matriz n X n la 

matriz unitaria. 


Para obtener 
1 
0 


0 


en la primera columna se requieren: n divisiones para producir el 1, 
y (n-—1)Xn operaciones para producir los 0. Esto significa que 
necesitaremos n + (n— 1) Xn =n? operaciones para conseguir la 
primera columna. Las operaciones que se requieren para obtener cada 
una de las columnas de una matriz unitaria n X n son también n?. 
Como hay n columnas, el número total de operaciones requeridas 
es n?, 


Habiendo encontrado la matriz inversa, tenemos que calcular A”! b, 
lo cual requiere otras n?* multiplicaciones. 
Así, el número total de operaciones que se requieren para obtener una 
solución es 

nn” +m, 
comparado con 


o 
3 3 


que exige el método de eliminación de Gauss. 


Número de | Eliminación de Gauss |Método de la inversa 
ecuaciones | Número de operaciones [Número de- operaciones 


36 

80 

1100 
=3 x 3,4x 10* 
=3 x 3,3x 108 


Cuando n es grande, el término n? es mucho más grande que el resto, 
de modo que el método de la inversa es, a grandes rasgos, tres veces 
más largo que el de eliminación de Gauss. Así, cuando hay que re- 
solver tres o más sistemas de ecuaciones simultáneas que tienen la 
misma A y diferentes b, es generalmente más rápido encontrar A”! 
primero, y utilizarla con todos los diferentes vectores b, en vez de 
calcular cada solución por el método de eliminación de Gauss. 


Ejercicio 1 


¿Cuántas multiplicaciones se requieren para efectuar la multiplica- 
ción de dos matrices n X n? |] 
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Ejercicio 1 
(3 minutos) 


MB 28.1.4 


28.14 Resumen 28.1.4 


En esta sección hemos medido la eficiencia, en términos del tiempo Resumen 
que requieren los cálculos, de tres métodos directos para resolver sis- ... 
temas de n ecuaciones (con n incógnitas). Para ello, hemos investiga- 

do el número de operaciones de multiplicar y dividir que requieren 

cada uno de ellos. Para valores grandes de n, hemos encontrado los 

siguientes resultados: 


Número 
aproximado de Utilidad cuando 
Operaciones n es grande 


El método 1,721 X En ningún caso es 
explícito aceptable. 
desarrollado en 

la Sección 28.1.1. 


Método de "| Ampliamente usado 

eliminación cuando no hay más 

de Gauss de tres sistemas 
con el mismo 
miembro de la 
izquierda. 


El método que Ampliamente usado 

utiliza la matriz cuando hay más de 

inversa 3 sistemas con el 
mismo miembro de 
la izquierda. 


Los métodos considerados en la Sección 28.1 son métodos básicos. 
Hay muchas formas especializadas de ellas aplicables a problemas 
particulares. 


Solución 28.1.3.1 


: 


Cada elemento del producto de dos matrices requiere n multiplica- 
ciones. 


Hay n? términos en la matriz producto. Por consiguiente, el número 
de multiplicaciones es n?. |] 


28.2 METODOS ITERATIVOS O INDIRECTOS 


28.2.0 Introducción 


Teóricamente, si se dispone de suficiente tiempo y se desprecian los 
errores de redondeo que pueden ocurrir en los cómputos, es posible 
resolver exactamente una ecuación matricial 


Ax =D, 


que tiene elementos exactos a; y b,, mediante alguno de los mé- 
todos directos descritos en la Sección 28.1. Por consiguiente, parece 
innecesario considerar cualquier otro método, como, por ejemplo, los 
métodos iterativos en los cuales conjeturamos una solución y la 
vamos refinando. Sin embargo, los métodos iterativos son importantes; 
los discutiremos por varias razones. Primero, con frecuencia sucede 
que, cuando hay que resolver un sistema grande de ecuaciones, un 
gran número de los elementos de la matriz de los coeficientes son 
iguales a 0; por ejemplo, 


xXx Xx 0 0 XX 0 
XK 0 Xx 0.00 
0%. x 0 0 0 
x 0.0 x 0 j 
<< 0 0 0. x-0 

Á 00 <xOO Ss 


donde cada X representa un número diferente de 0. Tales matrices 
se llaman matrices ralas (por razones obvias). Un método iterativo, 
que toma en cuenta los 0, puede, algunas veces, producir el resultado 
mucho más rápidamente que uno directo. Por supuesto, si la disposi.- 
ción de los O fuera conveniente, por ejemplo, si comenzáramos con 


N 


entonces la sustitución regresiva en el método de eliminación de 
Gauss sería un proceso directo. Lo importante es que el método ite- 
rativo no depende de una disposición conveniente. Segundo, el méto- 
do iterativo, si es convergente, mejora la exactitud de la solución 
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Solución 28.1.3.1 


28.2 


28.2.0 


Introducción 
* $ 


Definición 1 


* * 


aproximada, en cada paso. Así, se puede utilizar, una vez que se ha 
obtenido una primera aproximación mediante un método directo, 
para mejorar la exactitud de la solución. Puede haber otras ventajas 
que provienen de la economía de tiempo y de uso de espacio en un 
computador digital, pero esas ventajas necesitan un análisis más 
detallado. En términos de los cómputos automáticos, los métodos ite- 
rativos tienen frecuentemente la ventaja de que, en unos pocos pasos, 
repiten el mismo proceso relativamente sencillo, lo cual es ideal para 
una programación económica y eficiente. El método iterativo es tam- 
bién interesante desde un punto de vista puramente matemático: 
muestra cómo las ideas de la aplicación de intervalos de error numé- 
rico, introducidas en la Unidad 2, Errores y exactitud, se pueden ex- 
tender a la discusión de vectores. 


28.2.1 Algunos métodos iterativos 


Resumiremos el método iterativo que, para resolver ecuaciones defi- 
nidas en el conjunto de los números reales, tratamos en la Unidad 2, 
Errores y exactitud. Para resolver la ecuación 


x*-5x+3=0, 
reagrupamos los términos de modo que 
x = Hx? + 3), 
lo cual determina los elementos fijos o invariables del dominio de la 
función 
x—Hx? +3)  (xeR), : 
es decir, los elementos del dominio que son iguales a sus imágenes 


en el codominio. Una primera conjetura, xo, de la solución dio lugar 
a un nuevo número, 


So 4x3 $ 3); 


que esperábamos que estuviera más cerca de la verdadera solución, 
que llamaremos X. Los pasos siguientes formaron la sucesión cuyos 
elementos x, venían dados por 


x, = Hx7-1 +3), 
y que convergían a X, supuesto que la magnitud de un factor es- 
calar, determinado por la función 

x——Hx? + 3) (xeR), 
era menor que la unidad. 


Aquí, la ecuación definida en el conjunto de los números reales se ha 
sustituido por la ecuación matricial 


Ax — b=0. 


La solución, en vez de ser un número, es ahora un vector solución. 
La aplicación, que deduciremos de alguna reordenación de la ecuación 
matricial, tal como 


x = Gx + Hb, 


aplica ahora un espacio vectorial n-dimensional en un espacio vecto- 
rial n-dimensional. Buscamos un vector invariable, es decir, un vector 
que no sufre alteración en la aplicación. Buscaremos un algoritmo que 
_nos dé una sucesión de vectores que converge al vector solución. Por 
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28.2.1 


Discusión 
* 


analogía con nuestra experiencia previa en los métodos iterativos, 
también buscaremos algún número, correspondiente al factor escalar, 
que nos dé alguna sugerencia sobre la convergencia o divergencia de 
la sucesión que se investiga. Una considerable cantidad del análisis 
comprendido en esa investigación, está fuera de los objetivos de este 
curso, pero podemos obtener una sugerencia oportuna si examinamos 
un ejemplo elemental. Utilizaremos como ilustración un sistema de 
dos ecuaciones simultáneas con dos incógnitas, aun cuando, corrien- 
temente, el método solamente se utilizará en sistemas mucho mayores. 


Examinaremos el sistema de ecuaciones 
x, +4x,=6 
2x1 + x2=3 
Es fácil encontrar que las soluciones son 
Xq =2, a =L; 


esto es, el vector solución es 


() 


Por ahora vamos a posponer el enfoque matricial y sencillamente re- 
ordenaremos las ecuaciones en un par de formas. Ante todo, despe- 
jaremos x2, en la primera ecuación, en términos de x,, y en la se- 
gunda ecuación despejaremos x, en términos de x2, con lo cual 
obtendremos 


x, = 1,5 — 0,25x, 
x, = 2,5 — 0,5x2 


Ahora desarrollaremos las sucesiones x(?, x%? donde los elementos su- 
cesivos* están definidos por 


x* 1D = 1,5 — 0,25x( 
e = 2,5 e 0,5xP 


Supongamos que la primera conjetura del vector solución, 
xo 
o)' 


Obtenemos, entonces, 


es 


xQ) = 1,5 — 0,25 x 0 =1,5 
x4 = 25-05 xX0=2,5 


y la sucesión de los vectores solución estimados es, hasta el quinto 
término, 


N (5) ss [o e 
o)” l15)” loss)” 11,06)” 10,98 
*x En los procesos iterativos anteriores hemos utilizado el subíndice r, por ejemplo xr, 

para indicar la r-ésima estimación de la solución. Ahora, para evitar confusiones con 


los subíndices que ya se están usando y con el uso normal de los superíndices como 
exponentes, utilizaremos un superíndice entre paréntesis, por ejemplo x “?. 
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Intuitivamente, parece que esta sucesión converge el vector solución 
2 


As 

Ejercicio 1 

Use la reordenación 
x1 =6 — 4x, 


x2 =5-—2x,, 


0 
y comience con el vector ( , para obtener una sucesión de cinco vec- 


tores solución estimados, correspondientes al sistema de ecuaciones 
dado en el texto. ¿Parece ser una sucesión convergente? 7] 


Este último ejercicio indica, tal como podíamos esperar de nuestra 
experiencia previa, que algunas reordenaciones “trabajan” adecuada- 
mente, y otras no. Lo que quisiéramos tener es algunos medios de 
probar una reordenación, correspondiente al factor escalar discutido 
en la Unidad 2. Por ejemplo, querríamos probar 


Pe a + 0,5x1 E 2x2 
x2 = 2,5 — x1 +0,5x2 
que es más complicada que ninguna de las dos reordenaciones previas. 


Comenzaremos nuestra investigación de un “factor escalar” en las 
ecuaciones simultáneas examinando el proceso iterativo en términos 
más generales. Nuestras reordenaciones tienen la forma iterativa 


x"+*0D = Gx" + Hb 


donde G y H son matrices cuadradas. Por ejemplo, en la reordenación 
que nos sirvió para preparar la sucesión que parecía ser convergente, 
teníamos que 


XD = 2,5 — 0,5x 
qe =1,5 — 0,251 


Esto se puede escribir en forma matricial como 


e 0 —0,51 [xP 0 0,5|/6 
+1) = as a . 
xy —0,25 0 JixY 0,25 O MS 
Así, en este caso, 


0 —0,5 0 0,5 
(5= Y, H= ; 
0,25 0 0,25 0 


Ahora explicaremos lo que entendemos por convergencia de la suce- 
sión de los vectores solución estimados. Cuando estamos resolviendo 
una ecuación en R, podemos trabajar en términos de errores numéri- 
cos en un solo número y tratar de hacer este error lo más pequeño 
posible. Tratamos de encontrar una reordenación de la ecuación, y 
de ahí una función, que haga lo menor posible el intervalo de error 
en el que se encuentra la solución verdadera. Ahora estamos traba- 
jando en términos de vectores y parece natural examinar el compor- 
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Ejercicio 1 
(3 minutos) 


Discusión 
€ 


(continúa en la página 20) 


Solución 1 
(dd 1 o Es 
0 , 5 + =7 , 33 > =63 
No parece converger. E] 


(viene de la página 19) 


tamiento de un vector error. Definimos el vector error de la r-ésima 
iteración por 


e = ¿0 Y, 


donde x es la r-ésima aproximación al vector solución exacto X. 
Cuando la sucesión converge a X, el vector error se aproxima al vec- 
tor 0. Queremos que el vector error e? se haga “más pequeño”, en 
cierto sentido que aún no hemos definido, cuando r crece. 


El vector solución X debe satisfacer la ecuación 
X = GX +Hb, 


puesto que esta expresión no es más que una reordenación de la 
ecuación 


AX —b=0, 
del mismo modo que 
x = Hx? + 3) 
es una reordenación de la ecuación 
x?* —5x+3=0. 
De las dos ecuaciones 
x0+1D= Gx" + Hb 
X = GX + Hb, 


podemos obtener una relación entre los sucesivos vectores error. Te- 
nemos que 


xe+1 =-X=- Gx" =6X 
= G(s" — X), 
es decir, 
got = Ge”. 


En la próxima sección utilizaremos este resultado y estableceremos 
una medida para la convergencia. A primera vista puede parecer bas- 
tante fácil determinar si e“+!1 es “menor” que e“; por ejemplo, 
podemos decir que los elementos de e“+' deben ser todos más 
pequeños (en valor absoluto) que los correspondientes elementos de 
e“. Ahora bien, cuando reflexionamos sobre ello, no parece ser 
satisfactorio. ¿Qué podríamos decir si casi todos los elementos de 
e“+1) son menores que los correspondientes elementos de e(”, pero 
unos cuantos son ligeramente mayores? Mirando cada elemento 
individualmente no encontraremos un concepto aceptable: debemos 
buscar, como medida, un número único. 
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Solución 1 


Definición 1 
Rh k * 


Ejercicio 2 
Otras dos reordenaciones de nuestras ecuaciones originales 
x¡+4x,=6 
LA ES 
son 
(1) x, =6 — 4x, 
Xx2 =5 — 2x, 
(1) x, =3 +0,5x, — 2x2 
x2= 2,5 — x1 +0,5x2 


Escriba las correspondientes fórmulas de recurrencia en la forma 
matricial 


pnl) ES Gx" EL Hb, 
y calcule los cinco primeros vectores error en cada caso de la ecuación 


gr 1).= Ge", 


a 
comenzando por e = (>) ES] 


28.2.2 “Medición” de un vector 


Para estudiar la convergencia del método iterativo debemos tener 
alguna medida que nos permita comprobar si el vector error se está 
haciendo más y más pequeño en las sucesivas jteraciones. Podemos 
inventar una medida arbitraria cualquiera, y la convergencia puede 
muy bien depender de la medida escogida. Por tanto, debemos pensar 
un poco en la clase de medida que es razonable. 


En primer lugar, queremos que nuestra medida sea un solo número 
porque es fácil comparar números y decidir si la sucesión de tales 
números, obtenidos de los sucesivos vectores error, converge hacia 
el número que está asociado al vector error 0. Esto significa que, en 
el caso general, queremos definir una función que aplique R” a R. 


Comenzaremos por considerar algunas medidas que no son satisfac- 

torias, de modo que obtengamos una idea más clara de las condiciones 
, 

que debe verificar una “medida” conveniente. 


a; 0 
az 0 

Sean q4= > y 0= » |, el vector -0, 
A 0 


Definamos la función 


a—>a,  (aeR”, 


$ 


donde a, es una “medida” de a. 


Es evidente que esta medida no es satisfactoria porque, aunque es 
posible que la sucesión de los a,, obtenida de los sucesivos vectores 
error, converja a 0, esto no nos dice lo que está sucediendo a los otros 
elementos del vector error. 
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Ejercicio 2 
(4 minutos) 


28.2.2 


Tema principal 
hh * 


(continúa en la página 22) 


MB 28.2.1, 28.2.2 


Solución 28.2.1.2 Solución 28.2.1.2 
0 -4 O 
(1) a q y b 


Los cinco primeros vectores error son 


(> (Catholic) les 
BI” N-2al” 188)” V-16%) 16481 
Los vectores error se van haciendo mayores por cualquier ““me- 


dida” razonable, lo cual sugiere que no hay convergencia; esto 
está de acuerdo con nuestras conclusiones del ejercicio 1. 


05 -2 0,5 0 
A 
(0, a o Flo 05)? 


Los cinco primeros vectores error son 
al [0,50 — 268 2,259 — 28 3,1250 — 5,58 
(7) he - al de - e) eos + 3,1258)” 
7,0625a — 98 
Mi + 7,0625B)' 


No es fácil ver hacia dónde va esa sucesión. Por ejemplo, si toma- 
mos a =fB= 1, obtenemos los siguientes vectores error: 


(0 Moa El Bo aa 
1)" los)” 11,25)” -| 0,875)" 2,5625)" 
De nuevo tenemos una sucesión de vectores que no ofrece ningu- 
na esperanza. am 


(viene de la página 21) 


Consideremos la función 


a $ a (a e R”). 
Entonces, 
00, 


y queremos analizar si la sucesión de “medidas” que obtenemos de 
los sucesivos vectores error converge a cero. 


Pero tampoco esto es satisfactorio porque, por ejemplo, si n es par y 
los vectores error son 


22 


para todo r>N, donde N es algún entero positivo, entonces 
en — 0, pero el error de cada elemento de la solución estimada 
de nuestro sistema no es despreciable. 


La primera medida considerada no es aceptable porque no tiene en 
cuenta todos los elementos del vector error. La segunda tampoco lo es 
porque, aunque considera todos los elementos del vector, podemos 
llegar a un resultado falso si unos elementos son positivos y otros son 
negativos. 


Podemos mejorar la segunda medida. Para ello contamos con dos 
formas obvias: definimos las funciones 


n 1/2 
m:a-—| Y a!) (ae R”) 


i=1 


n 
m,:4— Y laj] (ae R”. 
i=1 
La medida definida por m, es la más obvia por dos razones; la me- 
dida es igual a la desviación estándar de a, de 0; además, en dos o 
tres dimensiones, esta medida se puede interpretar como la longitud 
o módulo del correspondiente vector geométrico. 


Analizaremos brevemente cada una de estas medidas para ilustrar 
su uso; nos limitaremos a la interpretación geométrica en dos o tres 
dimensiones. 


Examinemos nuevamente el ejemplo bidimensional que discutimos en 
la sección anterior y, en particular, la reordenación que nos llevó a la 
ecuación 


got = Gen 


0 = > 
G= : 
0,25 0 


Si escribimos 


donde 


entonces 

m2 — NP 
Si representamos el vector error mediante una flecha que vaya desde 
el origen de un conjunto de coordenadas cartesianas al punto cuyas 


coordenadas son (ef,ef), entonces m,(e“)) es la longitud de la 
flecha. Ahora, tenemos que 


A 
eg*o] —l-o25 0 Hep) 


esto es, 
ef +10 = —0,5e 


eg +0 = —0,25 ef. 
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. . . a ; 
Supongamos que el vector error inicial es g% = (7) entonces, pode- 


mos trazar una flecha desde el origen para representar el vector error. 
Los sucesivos vectores error están entonces, representados por 


Dirección 1 


En este ejemplo particular, el vector error tiene una de las dos direc- 
ciones y alterna entre ellas, de manera que, después de dos pasos 
en el proceso de iteración, el vector error está nuevamente señalando 
en la misma dirección, pero su longitud ha disminuido en un factor 8. 
Tomando un número suficiente de pasos, podemos hacer la longitud 
del vector error tan pequeña como se quiera, esto es, podemos con- 
seguir que la punta del vector esté tan cerca del origen como quera- 
mos. Esto sugiere que, con esta medida de un vector error, el método 


; 2d 4 A : a 
de iteración convergerá a una solución, cualquiera que sea el eros | 


de la conjetura inicial con la cual iniciamos la iteración. 


Ahora consideraremos la interpretación geométrica de la segunda 
medida, definida por 


ma) = Y al 


Usando la misma notación que antes, tenemos, en dos dimensiones, 
male”) = 11 + 109, 


que está representada como la suma de las dos longitudes que se han 
señalado con color en el diagrama. 


En nuestro ejemplo vemos que la suma de los módulos de los vectores 
error sucesivos se hace cada vez menor, de modo que una vez más 
tenemos un concepto intuitivo de que el vector error converge hacia 
el vector 0. 
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Dirección 2 


Dirección 1 


Como ya hemos conseguido una impresión intuitiva de lo que enten- 
demos por convergencia en los vectores, seremos ahora más precisos. 


Una medida que tenga características de las que hemos estado dis- 
cutiendo se llama norma. Una norma es una función, con ciertas pro- 
piedades especiales, que aplica los elementos de un espacio vectorial 
V en R. La imagen que, en esta función, tiene un vector v se denota 
por || || (y se lee “la norma de uv”). 


Se define que toda norma tiene las propiedades siguientes: 


(i) llell > 0, para todo pe V; 
(ii) lle] =0 si y solamente si y =0; 
(11) ll2, + v2ll < lle, ll + lp, ll, para todos p,, 1, € V; 
(iv)  [la«o,l| = la) Im, ll, donde «a es un número real cualquiera. 


Ahora definimos la convergencia en un espacio vectorial. Sea m una 
norma en el espacio vectorial y sea 


xD, ge plis 


una sucesión infinita de vectores. Se dice que esta sucesión tiene 
por límite a X si la sucesión de normas de los vectores error: 


ls —.X, 14? — X],... 


tiene por límite 0. Así, pues, hemos definido la convergencia en un 
espacio vectorial en términos de la noción de convergencia en R, que 
nos es familiar. No tenemos la intención de ampliar ahora el estudio 
de la normas. Un próximo paso obvio sería ver cuáles de nuestros re- 
sultados acerca de la convergencia en R se pueden extender a la con- 
vergencia en un espacio vectorial, pero esto lo dejaremos para un curso 
posterior. Ya hemos mostrado que podemos extender la idea de 
convergencia de una manera aparentemente satisfactoria, y que es 
suficiente para nosotros. En el próximo ejercicio tendrá que demos- 
trar un sencillo resultado acerca de la norma m2 introducida an- 
teriormente. 


En el trabajo que sigue supondremos que la norma es m2 , porque ha 
demostrado ser, en la práctica, más conveniente que m;. 


Ejercicio 1 
Considere la norma mz en R”, y demuestre que, si el límite de la 
sucesión 


A a es Xx 
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Definición 1 
* * * 


Notación 1 
+ + $* 


Ejercicio 1 
(5 minutos) 


entonces el límite de la sucesión 


DO U2) 3 k 
EI is io 
es X,, donde 
(k) 
X1 Xy 
(k) 
Xx X, 
na A "EE 
(k) 
a As 


Observe que este resultado nos dice que, si una sucesión de vectores 
es convergente, entonces las sucesiones de los correspondientes ele- 
mentos de los vectores mismos convergen a los correspondientes ele- 
mentos del vector límite. E 


Ejercicio 2 


Este ejercicio pretende dar una idea geométrica de m2 en dos di- 
mensiones. Más adelante nos referiremos a él. 


a e, 
Sea e= un vector error. 
e2 


(i) Trace la gráfica del conjunto solución de 
m(e) = le,] + lez] = 2 


(i.e. |x| + ly] =2 en coordenadas cartesianas). 
(SUGERENCIA: Considere, separadamente, los cuatro cuadran- 
tes del plano.) 
¿Qué implicación tiene este resultado sobre la flecha que parte 
del origen y que representa el vector error? 

(ii) ¿Qué región del plano representa el conjunto solución de 


le,| + lez] < 2? 


(e: lx + y] < 2). a 


Dada una sucesión de vectores error, tenemos ahora los medios para 
comprobar su convergencia por medio de la norma m2 en R”. Pero 
esto no es lo ideal porque tenemos que obtener la sucesión antes de 
poderla probar. Sabemos que la sucesión se obtiene de 


er+1) = Gel” 
si nuestra ecuación original 


Ax =b 


se reordena en la forma 
A) == Gx" El Hb. 


En otras palabras, la sucesión de vectores error está determinada por 
la reordenación; en particular está determinada por la matriz G. Por 
tanto, debemos determinar cuáles son las condiciones que podemos 
imponer a G (es decir, a la reordenación) de modo que la sucesión re- 
sultante de vectores error converja a 0. 


Ahora bien, ya hemos encontrado con anterioridad esta clase de pro- 
blema. En la Unidad 2, Errores y exactitud, resolvimos ecuaciones en 
R reordenándolas de modo que consiguiéramos una sucesión iterativa. 
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Ejercicio 2 
(4 minutos) 


Discusión 
* * 


Encontramos que existía un criterio asociado a la reordenación, el 
cual determinaba la convergencia: era una condición sobre el factor 
deescala. Examinemos nuevamente ese concepto. 


Supongamos que queremos resolver la ecuación 
f(x) =0, 

donde f es una función real y que tenemos una reordenación 
x = F(x) 

de esta ecuación, que nos lleva a la fórmula de recurrencia 
x0+1 = E(x), 

Definimos el factor de escala de la función F' como 


error estimado en la imagen de x 
. error en x 


es decir, 


error en x +1 
error en x 


Tomamos, por ejemplo, x(" como una conjetura de la solución de 
x = F(x), y calculamos el factor de escala, que suponemos que ha 
de variar muy poco en los alrededores de la solución. Si la magnitud 
del factor de escala fuera mayor que 1, entonces el intervalo de error 
que contiene la conjetura y la solución verdadera crecería, y si la 
magnitud del factor de escala fuera menor que 1, entonces el intervalo 
de error se contraería en cada aplicación de F y el proceso iterativo 
convergería hacia la solución. 


Veamos qué significa esto en nuestra discusión. Podemos definir 
aquí el factor de escala como 


ma(er +) _ norma del error estimado en x “+” 
ma(e”) norma del error estimado en x 


Observe que, por nuestra definición de norma, el factor de escala no 
puede ser negativo. Si este factor de escala es menor que 1 para todo 
r, entonces 


mae" +) < mae”), 
le. 
¡STD + 1eg7 0] + + [e +0] < 100] + le9] + --- + ]e(). 


Esta desigualdad se verifica para todo r, de modo que podemos en- 
contrar un número k, 0<kX <l, tal que 


n n 
bj jefr+2) < k de [ef”)] 
i=1 i=1 
n 


< k? AN lefr= 


i=1 


< ke Y le. 


i=1 
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Tema principal 
*h e +* 


(continúa en la página 29) 


Solución 1 
Puesto que el límite de x,x%,... es X, sabemos que la sucesión de 
números reales 
lx — Xi, lx — XI), ... 
converge a 0. 


Por consiguiente, usando nuestra definición de la norma mz, sabe- 
mos que el límite de la sucesión cuyo k-ésimo término es 


PP — Xd? — Xx — XK 
es 0. 


Ahora necesitamos un resultado obvio, pero que no ha sido demostra- 
do, que nos asegure que, si y y y son sucesiones cuyos términos son 
todos positivos, y si 


lím (u +.) =0 


entonces u y y convergen a 0. (Si lo cree conveniente repase sus co- 
nocimientos de la Unidad 7, Sucesiones y límites I, y demuestre este 
resultado.) 


Podemos dividir nuestra sucesión dada en n sucesiones, una de las 
cuales es 


a = Xp), 1x8 — X ls. + 
Puesto que todos los términos de la expresión 
PX iS A +12 


son positivos, podemos aplicarle el resultado que enunciamos más 
arriba. En vez de dos sucesiones, tenemos n sucesiones, una de las 
cuales se obtiene del primer término que aparece entre barras en la 
expresión dada anteriormente, esto es, 


1 2 
Ixg> E Xy, xP = Milos 
que converge, por consiguiente, a 0; esto es, la sucesión 
1d 2 
ea AA 


converge a X,. E 


Solución 2 


(i) En el primer cuadrante del plano, x > 0, y > 0 y el conjunto so- 
lución es 


xEy=2, 10. y =2=x: 
En el segundo cuadrante del plano, x < 0, y > 0, de modo que el 
conjunto solución es 

=x+y=2, 1e.y=2+x, 
y así sucesivamente. 


La gráfica del conjunto solución es, por consiguiente, el cuadrado 
que aparece en color en el diagrama. La punta de la flecha está 
sobre el cuadrado. 

(1i) La región que representa el conjunto solución es el interior del 
cuadrado coloreado. 
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Solución 1 


Solución 2 


(viene de la página 27) 


Se sigue que 


lím [3 $ 01) < [rs lím e) (E A) 
ande i=1 


r grandel ¡=1 
ya que > le] es constante. 
i=1 


Ahora, puesto que 0 <kXk < 1, lím k” = 0, de modo que 
r grande 
lím [5 ero] == () 
r grandel ¡=1 


Así, una vez más el factor de escala verifica nuestra orientación. Si el 
factor de escala es menor que 1, entonces la sucesión de normas de 
los vectores error converge hacia 0, es decir, la sucesión de los vec- 
tores error converge hacia el vector 0, y la sucesión de los vectores 
recurrentes x“? converge hacia la solución X. 


Tenemos, pues, un criterio para la convergencia, a saber, 


male" e 1) 
ma(e!”) 


Sabemos que 
ee) == Ge”, 
de modo que esa condición equivale a 


m,(Ge!”) 


== < 1, 
male!”) 


y parece que G verifica nuestra orientación sobre convergencia. La 
información contenida en G se relacionará, en alguna forma, con las 
normas de los sucesivos vectores error. 


Iustraremos el caso bidimensional. 
Si 
male" +1) < km,(e”), donde 0<k=<1, 
entonces las longitudes de los lados de los cuadrados disminuyen en 
un factor k en cada paso y, en consecuencia, la flecha que representa 


el vector error y que está contenida en el cuadrado apropiado se aproxi- 
ma al vector cero, es decir, la sucesión de vectores es convergente. 
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(Véase el ejercicio 2.) 


e; 


Examinemos ahora lo que significa la condición del factor de escala, 
en términos de los elementos de G. Realizaremos el análisis solamente 
en dos dimensiones, pero los argumentos empleados se pueden genera- 
lizar fácilmente. 


Supongamos que 


E ls E, 
821 822 


Los vectores error sucesivos son, entonces, 


[eo] = A nl 
e 821 £2211e7 
de modo que 
+1) 


fr == (r) (r) 
ei = Buitre + 201202 


(+1) = (r) (r) 
ez = 821801 + 822€2- 


maeotD) = 10901 + 10579 
= 181100 + 81228 + lez," + 82209) 
< lg81, 80) + lg] + lg2,8 + 1222891, 


de acuerdo con la desigualdad triangular (véase la Unidad 22, Algebra 
lineal I, Sección 22.1.3). 


Puesto que, para cualesquiera números reales «, $, 
la x Bl =]la] x 18, 
tenemos que 
(A +19 < lg 110 +18, 2119 + 18211180 +l82211e2), 
finalmente, 
¡STD + 180 < (8111 +l8211)1e71 + (28121 + 8221) 189. 
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Por consiguiente, si podemos encontrar un número k en el intervalo 
0<k<l1 tal que 


Igual + lgz1l <k 


y que, además, 


lg12l + l8221 < k, 


se verifica la condición del factor de escala, y la sucesión de vectores 
error converge hacia el vector 0. Así, podemos probar la matriz G su- 
mando los módulos de los elementos de cada columna. Si la mayor 
suma obtenida en cualquier columna es menor que 1, la sucesión de 
vectores error es convergente. Por ejemplo, si 


po el 
—Noj2 0,1)” 


entonces tenemos que 


gil +1l8211=0,5 y  lgrl + lgz21 = 0,6 


Si tomamos k =0,6 (o cualquier número entre 0,6 y 1), entonces se 
verifica la condición del factor de escala, de modo que la sucesión de 
recurrencia producida por la matriz G es convergente. 


(Este número k, que representa la mayor de las sumas de los módulos 
de las columnas de una matriz, puede, a su vez, tomarse como una 
norma de una matriz. En otras palabras, es un método para “medir” 
una matriz o para asociar un número a una matriz.) 


Ejercicio 3 
Pruebe en cada una de las siguientes matrices si la sucesión de vec- 
tores error obtenida de la ecuación 

gr+no = Ge" 


converge realmente hacia el vector 0. Si es así, encuentre un valor 
conveniente para la constante k. 


] 0,3 0,75 
Gi) G= 
0,6 0,2 
A 0,7 0,05 
(11.6 = 
—0,5 0,1 
0,4 0,3 
(110) 6 = 
0,6 0,7 
0,33 —0,72 
(in G= 
—0,66 —0,27 s 


Ejercicio 4 
Reordene las siguientes ecuaciones de tal manera que el método itera- 
tivo resultante sea convergente. Después, resuelva las ecuaciones por 
iteración 

Sx, + 3x, = 7 

2x, — 4x, = 3. li] 
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Ejercicio 3 
(3 minutos) 


Ejercicio 4 
(3 minutos) 


(continúa en la página 32) 


Solución 3 


(i) La sucesión es convergente. Cualquier k tal que 0,95 <k < 1 es 
conveniente. 
(11) La sucesión no es necesariamente convergente. 
(iii) La sucesión no es necesariamente convergente. 
(iv) La sucesión es convergente. Cualquier k tal que 0,99 < 
conveniente. 


A 


< l es 


Observe que solamente hemos mostrado que 

si se verifica la condición => existe la convergencia 
pero no hemos demostrado que 

si no se verifica la condición = existe la convergencia. 


Es posible que aún en este último caso exista la convergencia, pero 
sería más difícil probarlo. Así, sería conveniente que examinara la 
sucesión generada por las G en los casos (ii) y (iii) comenzando con 


() . 


Solución 4 
Una reordenación es 
x1 = 1,4 — 0,6x2 
x2 = —0,75 +0,5x 1 
Con esta reordenación, 
( le 
G= 
0,5 0 
y la fórmula de recurrencia toma la forma 
yrD ( ES a 1,4 ) 
0,5 0 —0,75 


Comenzando con cualquier x(', llegaremos, eventualmente, a estar 


37 

26 1,42 
tan cerca como queramos de la verdadera solución i = 60d) 

26 Ñ 


(viene de la página 31) 


La prueba del factor de escala se puede aplicar a una sucesión de vec- 
tores error en un espacio vectorial de cualquier dimensión finita. 
Sumamos los módulos de los elementos de cada columna de una 
matriz y si todas las sumas son menores que o iguales a k, que es un 
número positivo < 1, entonces hemos garantizado la convergencia. 
Por ejemplo, para el caso de tres dimensiones, tenemos otra interpre- 
tación geométrica. Usando la norma le,| + lez| + lez], encerra- 
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Solución 3 


Solución 4 


Discusión 
* * 


mos la flecha que representa el vector error dentro de un octaedro que, 
en cada paso, se hace menor si el método converge. 


Para más de tres dimensiones carecemos de representaciones geo- 
métricas pero el álgebra es la misma. 


Por último, veamos un punto que mencionamos en la introducción de 
esta sección. Dijimos que el método iterativo era útil cuando la matriz 
de los coeficientes es una matriz rala. En efecto, cuando en cada fila 
hay solamente unos pocos elementos diferentes de 0, podemos fácil- 
mente aislar uno de los elementos del vector solución y expresarlo en 
términos de relativamente pocos elementos del miembro de la derecha. 


(Quizá sean necesarias algunas manipulaciones para obtener tal ecua- 
ción para cada elemento.) Por ejemplo, la reordenación 

Xi =.0,2%, 0,5xs + 1 

X2 = 0,4x, +0;1%3 + 2 

X3 = 0,2x, E 0,3X4 + 3 

Xa = 0,5x3 +0,2x5 +4 

x5 =0,4x, +0,3x2 + 5 
tiene 3 multiplicaciones menos en el miembro de la derecha de cada 


ecuación, que las que debería tener. Sucede así porque la matriz de 
los coeficientes escrita en la forma 


1 =0Z2 0 O  =05 
—0,4 il, 0,1 0 0 

oO  —0;2 1 0.3. 0 

0 0: + =0/5 LL — =052 
07403 0 0 1 


es relativamente rala. Se notaría aún más la notable disminución del 
número de multiplicaciones requeridas si, por ejemplo, solamente hu- 
biera cinco elementos diferentes de O en cada fila de una matriz 
100 Xx 100. Sin embargo, la decisión final sobre la eficiencia, com- 
parada con la de un método directo, es más difícil porque tenemos 
que decidir también sobre la rapidez de convergencia del proceso ite- 
rativo, supuesto que converge; esto es, determinar cuántos pasos se 
requieren para llegar a la exactitud requerida. 


33 


MB 28.2.2 


Ejercicio 5 


El siguiente es un sencillo programa rutinario, llamado $ ITER, que 
se puede usar para obtener soluciones de ecuaciones simultáneas 
por iteración. Usted debe ser capaz de comprender cómo se desarrolla 


este programa, por simple inspección. 


95 
100 
105 
110 
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PROGRAMA 
PRINT “X = AX + B, A ES UNA MATRIZ N*N” 
PRINT N =” 
INPUT N 
IF N < = 10 THEN 20 


PRINT “N ES DEMASIADO GRANDE; DE UN NUMERO 
DE 1 A 10” 


GO TO 10 


PRINT “A =>” 
FOR I=1 TO N 
FOR J=1 TO N 
INPUT A(I, J) 
NEXT J 

NEXT 1 


PRINT “B =>” 
FOR I=1 TO N 
INPUT B(I) 
NEXT I 


PRINT “X INICIAL =>” 
FOR I=1TO N 

INPUT X(I) 

NEXT 1 


PRINT “¿QUE SIGUE? RESPONDA CON —1 SI SE 
NECESITAN MAS ITERACIONES” 


PRINT “+1 SI SE VA A INSERTAR UNA NUEVA MATRIZ 


Y 0 PARA TERMINAR” 

INPUT X 

IF X = —1 THEN 110 

IF X = +1 THEN 5 

IF X =0 THEN 300 

PRINT “LA RESPUESTA DEBE SER -—1, +1 6 0” 
GO TO 80 

PRINT “¿CUANTAS ITERACIONES?” 
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Ejercicio 5 


COMENTARIO 


Prepara el trabajo. 


Pregunta por las 
dimensiones de la 
matriz, que debe ser 
menor que 10. 


Da entrada a la matriz 
A, fila por fila, un 
elemento cada vez. 


Da entrada a B. 


Pide la estimación 
inicial de X. 
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PROGRAMA COMENTARIO 
115 INPUT M Asegura que no se 
120 IF M< 21 THEN 135 necesitan más de 20 
iteraciones. 


125 PRINT “DE UN NUMERO DE 1 A 20” 
130 GO TO 115 

135 FOR K=1 TO M 

140 FOR I=1 TO N 

141 YM) =B(D) 

142 NEXT 1 

145 FOR I=1 TO N 

150 FOR J=1 TO N 


155 IF A(L, J) = 0 THEN 165 Efectúa las iteraciones. 
160 LET Y(D) = YM) + A(L, J) *X(J) 

165 NEXT J 

170 NEXT 1 


175 FOR I=1 TO N 

176 XM) = YM) 

177 NEXT I 

180 NEXT K 

185 PRINT “X =>” 

19 FOR I=1 TO N 

192 PRINT X(D, Imprime los resultados. 
194 NEXT 1 

196 GO TO 75 

300 END 


(i) Tal como está presentado, este programa no puede decirnos si la 
matriz A producirá una iteración convergente. Escriba un conjunto 
de instrucciones que modifiquen el programa de tal manera que 


n 
nos diga el máximo valor de la] para ¡=1,...,n, y no efec- 
. a 


túe la iteración a no ser que este valor sea menor que 1. 
(11) Utilice el programa para resolver el siguiente sistema 5 X 5 con 
dos cifras decimales. 


x, — 0,2x2 —0,5x5 = 1 
—0,4x1 + 1,1x2 — 0,1x3 =2 
—0,2x2 + A 0,3x4 =3 


—0,5x3 + 1,2x4—0,2x5 = 4 
—0,4x1 = 0,3x2 +3 Xs= S: |] 
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Solución 5 


() 


(11) 


Un conjunto de instrucciones adecuadas sería: 


46 GOTO 200 
200 LET P=0 
205 FOR J=1 TO N 
210 LET Z=0 


215 FOR I=1 TO N 

220 LET Z =Z + ABS(A(I, J) 
225 NEXT I 

230 IF Z <P THEN 240 
235 EETP =Z 

240 NEXT J 


245 PRINT “NORMA DE LA MATRIZ =”; P 


250 IF P < 1 THEN 50 


255 PRINT “NORMA DE LA MATRIZ DEMASIADO GRANDE” 


260 GOTO 76 


Una reordenación conveniente es 
0 0,2 0 0 0,5 
0,4 —0,1 0,1 0 0) 


x=| 0 0,2 0 0,3 0 lx+ 


0 0 0,5 —0,2- 0,2 
0,4 03 0 0 0 
La solución, con dos cifras decimales, es 
6,44 
4,72 
x=| 6,16 
7,40 
8,99 


28.2.3 Resumen 


Para resolver el sistema de ecuaciones representado en forma matri- 


cial por la ecuación 


Ax = b, 


la reordenamos en la forma 
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x = Gx + Hb. 


=— 


HH bh uu nn 
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Solución 5 


28.2.3 


Resumen 
* * 
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Esto nos lleva a una posible sucesión recurrente definida por 
A+D=G" + Hb  (r=1,2..,) 


que será convergente si la suma de los módulos de los elementos de 
cada columna de G es menor que 1. Comprobamos este resultado 
examinando el comportamiento de los vectores error 


¿€ =x0-X (P=1,2. 2) 
que verifican 
go+D = Ge” (r ==, Le: . =d 


donde X es el vector solución exacto. 


28.3 SISTEMAS DE ECUACIONES DE MAL 28.3 
COMPORTAMIENTO 


En esta sección examinaremos una forma particularmente peligrosa Tema principal 
que ocurre, algunas veces, cuando resolvemos ecuaciones simultáneas Es 
basadas en datos inexactos o cuando introducimos errores de redon- 

deo durante la solución. De hecho, la acumulación del error puede 

hacer virtualmente inútiles los resultados de la solución de un sistema 

de ecuaciones simultáneas. Cuando un pequeño cambio en los datos 

producen un efecto grande en la solución de un sistema de ecuaciones, 

decimos que el sistema es de mal comportamiento. Este término no 

tiene una definición precisa. “Pequeños cambios en los datos” que 

producen “grandes efectos en el resultado” no es más que una ma- 

nera de describir el fenómeno del mal comportamiento; los adjetivos 

usados indican la imprecisión del concepto. En el próximo ejemplo 

daremos una ilustración numérica del mal comportamiento. 


Ejemplo 1 Ejemplo 1 


En la reconstrucción de un crimen se encontró que una bala había 
perforado en A y Buna ventana de doble vidrio. Las distancias, medi- 
das desde el centro de cada perforación, se indican en el diagrama, 
y están aproximadas al centímetro más cercano. Otra evidencia de- 
mostró que el arma se encontraba en el nivel CD cuando fue dispara- 
da. ¿Con cuánta exactitud se puede determinar la posición desde 
donde se disparó el arma? 


Xx 
A Y e 
20cm 
D 
L— 10 cm —» X; A 


Antes de continuar con la solución de este ejemplo, le será útil supo- 
ner que los 10cm han sido medidos exactamente, y trazar cuidadosa- 
mente, sobre una hoja de papel milimetrado, las dos posiciones extre- 
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mas de las posibles trayectorias de la bala, supuesto que éstas son 
líneas rectas. 


Solución del ejemplo 1 


Si hacemos coincidir con el diagrama unos ejes coordenados x, Y Xz, 
estaremos realmente buscando la intersección de la recta AB con 
el eje x,, CD. 


Hojas de vidrio 


Una posible 
trayectoria de la bala 


e D 


La ecuación de la recta AB es 
X2 = mx; + C, 


donde (teóricamente) podemos encontrar m y c a partir de los datos. 
La distancia xc, desde C hasta el punto de donde salió la bala, es 
el valor de x, cuando x2 = 0, es decir, 
e E 
Ji 
Si suponemos que los datos son exactos, tendremos que AB pasa por 
los puntos (0, 20) y (10, 18), de modo que 


20 =0 
18 = 10m + c, 
20 
de donde m = —0,2. Por tanto, xa = cm = 100 cm = 1 m. 


0,2 


Ahora bien, los datos no son exactos: c podría ser 19,5 y, entonces, m 
podría ser tan pequeño como nos dijera la ecuación 
18,5 = 10,5m + 19,5 

1.e. 

m= es 

10,5" 

de modo que x¿« podría ser tan grande como 19,5 X 105cm = 
204,75 cm. 


Esto es, el valor aproximado de xc (1m) puede contener un error de 
más de un metro. En verdad, todo lo que podemos decir es que 


xa E [64,92, 204,75] cm 


En otras palabras, la distancia desde donde fue disparada el arma 
(medida desde la ventana) podría ser cualquier valor desde 0,5 m 
hasta 2 m. 


Este resultado no es de gran precisión si consideramos la aparente 
precisión de las mediciones originales. la] 


En términos geométricos, en el último ejemplo estuvimos intentando 
encontrar la intersección de un par de rectas casi paralelas (el eje x, 
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(Véase RB3) 


Discusión 
* * 


era una de las rectas). Si hay un pequeño error en los datos originales, 
puede producirse un error considerable en la solución. Un modo de 
imaginarnos lo que está sucediendo consiste en considerar dos lin- 
ternas cuyos rayos de luz se cruzan. 


La solución exacta, si los datos son exactos, está representada por el 
punto P. Si los datos son inexactos, hecho que se representa por los 
rayos divergentes que salen de las linternas, entonces las rectas 
podrían ser cualesquiera de las que forman los haces luminosos. En 
consecuencia, el punto que representa la verdadera solución podría 
ser cualquiera de los que forman “el área de intersección” que aparece 
coloreada en el diagrama. 


Intuitivamente, cuanto más tiendan a ser paralelos los ejes de las 
linternas, tanto más alargada será el área sombreada y, consecuente- 
mente, tanto más lejos de la verdadera solución se puede encontrar 
la solución aproximada. 


Examinemos un caso más general, en dos dimensiones, en el ejercicio 
siguiente. 


Ejercicio 1 


Resuelva, por el método de eliminación de Gauss, las siguientes ecua- 
ciones simultáneas: 


xi + x=1l 
xy + (1 + ex, =2, 


donde «e es algún número real. En las fórmulas que obtenga para xi 
y x2 sustituya para 


(i) e = 0,01, 0,02 
(ii) e = 2,01, 2,04 


respectivamente. Compare los cambios en porcentaje del coeficiente 
de x2 en la segunda ecuación, en los casos (i) y (ii), con los cambios 
en porcentaje de la correspondiente solución de xz. E 


Este último ejercicio es más bien hipotético. Sin embargo, los errores 
de redondeo que ocurren en los datos iniciales de un problema prác- 
tico y en el proceso de solución mismo, no son hipotéticos. Estos pe- 
queños errores de redondeo pueden conducir fácilmente, en ecuacio- 
nes de mal comportamiento, a resultados muy inexactos. Un modo 
sencillo de probar el mal comportamiento consiste, precisamente, en 
proceder como lo hicimos en el último ejercicio, esto es, efectuar un 
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Ejercicio 1 
(4 minutos) 


Discusión 
* $ 


pequeño cambio en algunos coeficientes para ver qué efectos se pro- 
ducen en la solución, pero esto es difícil de hacer cuando se trata de 
conjuntos de ecuaciones mayores que el considerado aquí, y parti- 
cularmente cuando una solución parece aceptable desde el punto de 
vista físico. Hay otros métodos que nos dan una indicación práctica 
cuando se presenta el mal comportamiento pero la mayor parte de 
ellos comprenden muchas y complicadas técnicas. Describiremos 
solamente uno de tales métodos. 


Teóricamente, el mal comportamiento se puede describir de la si- 
guiente manera. Dado un sistema de ecuaciones simultáneas en 
forma matricial 


Ax = d, 


el sistema es de mal comportamiento si las n columnas de la matriz 
son casi linealmente dependientes (nótese nuevamente la vaguedad 
de los términos) o, en otras palabras, si la matriz es casi singular. 
Esto significa que todos los elementos de A están cerca de los corres- 
pondientes elementos de alguna matriz singular. En el último ejer- 
cicio, por ejemplo, A es 


ñ | 
l 1+€ 
Evidentemente, si e fuera pequeño (ya vimos que esa es la condi- 
ción que convierte las ecuaciones en ecuaciones de mal comporta- 


miento), un cambio pequeño, es decir, un cambio de —e en 022, 
convertiría la matriz en una matriz singular, 


l, 
¡| 
en la cual las dos columnas son obviamente linealmente dependientes. 


Ejercicio 2 


Determine la matriz inversa A”! de 


Al 
7 ; 
l 1+e 


¿Cuáles son los elementos de A7! si e = 0,01? ll 


Este último ejercicio muestra una característica de la inversa de la 
matriz de los coeficientes de un sistema de ecuaciones de mal com- 
portamiento: sus elementos son relativamente grandes cuando se 
comparan con los elementos de la matriz original. Esto nos lleva a una 
observación final acerca de los sistemas de ecuaciones que pueden 
ser de mal comportamiento. Un medio recomendado y útil para probar 
la exactitud de una solución estimada x, de 


Áx =b 
es multiplicar x,, a la izquierda, por A y obtener un vector b,. Esto es, 
Ax, = Db,. 


Si los elementos de b, difieren muy poco de los elementos de b, en- 
tonces podremos suponer, normalmente, que la solución es aceptable- 
mente exacta. (Recuerde que casi siempre habrá algunas inexactitu- 
des en los cómputos, cualquiera que sea su extensión, provenientes 
de los errores de redondeo.) Si el sistema de ecuaciones es de mal 
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Ejercicio 2 
(4 minutos) 


Discusión 
* $ 
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comportamiento, sin embargo, las soluciones pueden contener mayores 
errores aún. En efecto, si sustraemos las dos ecuaciones anteriores 
obtenemos 


A(x, — Xx) = b, — b, 


o, usando la notación del vector error, 


donde 


e=x-x y E=b,-b 
Entonces tendremos 
e= AE. 


Si los elementos de A”! son grandes, es intuitivamente obvio que e 
puede ser grande aun cuando E sea pequeño. Así, siempre que sos- 
pechemos que hay mal comportamiento, vale la pena comparar los 
elementos de A7* con los elementos de A. 


Ejercicio 3 Ejercicio 3 


(4 minutos) 
Si 
Lo A 
A=l, a 
1. 10% 


encontrar un E con una norma (suma de los módulos de los elementos) 


< 0,01 que conduzca a un e cuya norma sea > l. Eu 
Finalmente, existe aún el problema de qué hacer acerca del mal Discusión 
comportamiento cuando existe y lo hemos identificado. Puede ser «$ 


_ conveniente un replanteamiento del problema en términos de otras 
variedades, como puede hacerse quizá en los conjuntos de ecuaciones 
simultáneas de mal comportamiento que aparecen en los problemas 
de estructuras en ingeniería, de modo que se llegue a un conjunto de 
ecuaciones de buen comportamiento. De lo contrario, podemos utilizar 
operaciones de doble precisión en los cómputos, esto es, trabajar con 
el doble número de dígitos, en un intento por obtener mejores resul- 
tados, pero si las ecuaciones son de extremado mal comportamiento, 
tampoco conseguiremos resultados aceptables. Entonces, sólo queda 
un consejo: olvídese del problema o acepte el resultado, cualquiera 
que sea su exactitud, aunque sea muy mala. 


Resumen Resumen 
* * 


Cuando la matriz de los coeficientes es casi singular, esto es, cuando 
un cambio pequeño en algunos elementos de la matriz produce una 
matriz singular, la ecuación matricial 


Ax =b 


es de mal comportamiento. Esto generalmente significa que la solución 
es altamente dudosa, en particular si los datos originales, de donde 
se formó la matriz, eran inexactos. 
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Solución 1 
Se x=1 
x1 + (1 + e)x, = 2 
Aplicando el método de eliminación de Gauss obtenemos 
x>+x=1l 
ex =1, 


de donde, 
1 
E le dal 0), 
€ 


y, por sustitución regresiva, 


1 
41 


(1) 


0,01 0,02 2,01 2,04 


(—99, 100)| (—49, 50) ((0,50, 0,50)| (0,51, 0,49) 


El cambio en porcentaje del coeficiente de x2 es 1% con una cifra 
decimal en cada caso. 


Solución con dos 
cifras decimales 


En (i) la solución x2 varía en un 50%; en (ii) en un 2%. Esto indica 
que “cuando e es pequeño las ecuaciones son de mal comporta- 
miento” significa, en este caso, “pequeño comparado con la unidad”. 


li 
Solución 2 
Es E 
q € € 
2. A 
€ € 
101 —100 
los 100 E 
Solución 3 
Un ejemplo es 
—0,004 
> cdo) 
Entonces, m2(E) = 0,008 < 0,01 y 
101 —100)| [—0,004 — 0,804 
ion 1o)l 0004) = | 0,8 
con male) = 1,604. || 
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Solución 1 


Solución 2 


Solución 3 


28.4 CONCLUSION 


En este texto hemos examinado un problema práctico: cómo resolver 
los sistemas de ecuaciones lineales simultáneas, particularmente los 
sistemas grandes. 

En la Sección 28.3 prestamos atención al fenómeno del mal compor- 
tamiento que puede aparecer en dichos sistemas y que pueden hacer 
muy difícil, y aun imposible, hallar una solución de las ecuaciones 
que sea aceptable o útil. 

Hemos demostrado los métodos básicos, directos e indirectos, de so- 
lución. Existen muchas variedades particulares de estos métodos. 
No obstante, lo único que podemos hacer aquí es examinar los méto- 
dos básicos y discutir su eficiencia y exactitud. 


Posdata 


“Esto es llevar las cosas muy lejos; yo termino aquí”. 


Harry B. Smith 
We Draw the Line at That 
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Conclusión 
* * 


M 100 — UNIDADES DEL CURSO BASICO DE MATEMATICAS 


00315004 0N-A 


Funciones 

Errores y exactitud 
Operaciones y morfismos 
Diferencias finitas 

NO HAY TEXTO 
Desiguales 

Sucesiones y límites 1 
Computación Il 

Integración I 

NO HAY TEXTO 

Lógica 1 — Algebra de Boole 
Diferenciación I 
Integración II 

Sucesiones y límites II 
Diferenciación HI 
Probabilidad y estadística 1 
Lógica II — Prueba 
Probabilidad y estadística II 
Relaciones 

Computación II 
Probabilidad y estadística II 
Algebra lineal 1 

Algebra lineal II 
Ecuaciones diferenciales 1 
NO HAY TEXTO 

Algebra lineal II 

Números complejos 1 
Algebra lineal IV 

Números complejos II 
Grupos I 

Ecuaciones diferenciales II 
NO HAY TEXTO 

Grupos II 

Sistemas numéricos 
Topología 

Estructuras matemáticas 
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